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Abstract 

In this note we study a minimization problem for a vector of measures subject to a prescribed interaction 
matrix in the presence of external potentials. The conductors are allowed to have zero distance from each 
other but the external potentials satisfy a growth condition near the common points. 

We then specialize the setting to a specific problem on the real line which arises in the study of certain 
biorthogonal polynomials (studied elsewhere) and we prove that the equilibrium measures solve a pseudo- 
algebraic curve under the assumption that the potentials are real analytic. In particular the supports of 
the equilibrium measures are shown to consist of a finite union of compact intervals. 



1 Introduction 

In this short paper we consider a vector potential problem of relevance in the study of the asymptotic 
behavior of multiple orthogonal polynomials for the so-called Nikishin systems [I]. The original 
problem was introduced in [2] (without external fields) and further questions has been addressed in 
[31 m [5l [H |7] . The main motivation of our interest for this problem arises in a recently introduced 
set of biorthogonal polynomials [8J. These polynomials are related on one side to the spectral theory 
of the "cubic string" and the DeGasperis-Procesi peakon solutions of the homonymous nonlinear 
differential equation [9J; on the other end they are related to a two-matrix model [lOj with a measure 
of the form ^ ^ , 

where the Mj's are positive definite Hermitian matrices of size NxN, a, (3 are some positive densities 
on M+ and the expressions a(Mi),/3(M2) stand for the product of those densities on the spectra of 

The relation between the relevant biorthogonal polynomials and the above-mentioned matrix 
model is on the identical logical footing as the relation between ordinary orthogonal polynomials and 
the Hermitian random matrix model [Ifj. 
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In [8] a Riemann-Hilbert formulation (similar to the formulation of multiple-orthogonal polyno- 
mials as explained in [12J but adapted to the peculiarities of the model) was derived and in (TU] the 
correlation functions of the spectra of the two matrices were completely characterized in terms of 
the matrix-solution of that Riemann-Hilbert problem. 

In [13J the analysis of the strong asymptotics with respect to varying weight (following [14j) will 
be carried out. A pre-requisite of that analysis is the existence and regularity of the solution of a 
suitable potential problem, namely the one which we explain in the second part of the paper. 

In fact, the present paper is addressing a wider class of potential problems that will be necessary 
for the study of the spectral statistics in the limit of large sizes of the multi-matrix model 

corresponding to a chain of positive-definite Hermitian matrices Mj with densities aj as above. 

In Section [2] we consider the problem as a vector potential problem in the complex plane with a 
prescribed interaction matrix. Under a suitable growth condition for the external potentials Vj{z) 
near the overlap region of the conductors (in particular the common points on the boundaries) it is 
shown that the minimizing vector of equilibrium measures has supports for the components separated 
by positive distances. 

In Section[4]we specialize the setting to the situation in which the conductors = (— 1)-'~^[0, oo) 
(so that they have the origin in common), with an interaction matrix of Nikishin type as in [T2]. We 
prove that the minimizing measure is regular and supported in the interior of the condensers (under 
our assumption of growth of the potentials). 

This result allows to proceed in Section [5] with a manipulation of algebraic nature involving the 
Euler-Lagrange equations for the resolvents (Cauchy transforms) Wj{x) of the equilibrium measures. 
It is shown that certain auxiliary quantities Zj that depend linearly on the resolvents and the potentials 
(see (5-4) for the precise formula) enter a pseudo-algebraic equation of the form 

+ C2{x)z^-^ + ... + Cr+i{x) = (1-3) 

where the functions Cj{x) are analytic functions with the same singularities as the derivative of 
the potentials V^'(x) in the common neighborhood of the real axis where all the potentials are real 
analytic. In particular the coefficients Cj{x) do not have jumps on the real axis and the various 
branches of eq. 1-3 are precisely the Zj{x) defined above. For example, if the derivative potentials 
are rational functions, then so are the coefficients of (1-3). This immediately implies that the 
branchpoints of (1-3) on the real axis (i.e. the zeroes of the discriminant) are nowhere dense and 
hence a priori the supports of the measures must consist of a finite union of intervals (since they 
must be compact as shown in Sect. [2] in the general setting). 
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The role of the pseudo-algebraic curve (1-3) is exactly the same as the well-known pseudo- 



hyperelliptic curve that appears in the one-matrix model [TBI EZ]: in the context of the study of 
asymptotic properties of multiple orthogonal polynomials it has been pointed out since the funda- 
mental work |[l5| that the Cauchy transforms of the extremal measures solve an algebraic equation. 

We also mention the recent work [18J, in which the limiting behaviour of Hermitian random 
matrices with external source is investigated and the presented asymptotic analysis relies on a set 
of conditions which are shown to be equivalent to the existence of a particular algebraic curve. The 
methods used in that paper to prove the existence for some special cases are very similar to our 
approach. 

As it was pointed out by one of the referees, examples of of algebraic curves for special external 
fields were also obtained in the recent papers [71 119]. 

1.1 Connection to a Riemann— Hilbert problem 

The principal motivation to the present paper is the application to the study of biorthogonal (multiply 



orthogonal) polynomials that arise in the study of the model hinted at by eq. (1-1). In |8l HOj we 
introduced the biorthogonal polynomials 

r e-NiVr{x)+V2iy)) 

/ Pn{x)qmiy) , dxdy = Cj,6mn , Pnix) = x"- + . . . , qn{y) = y"' + . ■ . . (1-4) 

JrI x + y 

In [T0| it was shown how a natural vector potential problem (for two measures) arises in that context 



and leads to a three-sheeted spectral curve of the form (1-3). Such problem enters in a natural way 
in the normalization of the 3x3 Riemann-Hilbert problem considered in [8] characterizing those 
polynomials (and some accessory ones) in the limit ^ oo, n — > oo, ^ — > T > 0. The notation 
Vi, V2 is meant here to reflect the notation that will be used in Sect. [Sjand Sect[6](up to a reflection 
V2{y) ^ V2{-y), as explained in U^)- 

In perspective the more general situation with several measures considered in Sect. |4]and Sect, 
[sjwill be associated to the polynomials appearing in the study of the random-matrix chain ( |l-2[ ) and 
biorthogonal polynomials for pairings of the form 



/ Pn{xi)qm{xK) ^^_^^ -\\dXj = c\5mn (1-5) 

nj=i (^i + ^i+i) i=i 

Pn{x) = + . . . , qn{y) = + . . . (1-6) 
The details are to appear in forthcoming publications [13] . 
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2 The vector potential problem 



In this section we consider the vector potential problem which is a slightly generalized form of the 
weighted energy problem of signed measures ([20], Chapter VIII; [4j, Chapter 5). 

Let A = {aij)^j^i be an R x R real symmetric positive definite matrix (in particular it has 
positive diagonal entries), referred to as the interaction matrix, containing the information on 
the total charges of the measures and their pair interaction coefficients. Suppose Si, S2, . . . , S/j is 
a collection of non-empty, not necessarily disjoint closed subsets of C such that T^k n has zero 
logarithmic capacity whenever < 0. Define the functions hf^-.C —>■ (— cxd,cxd] for each T^j. to be 

/i,(z):=ln— , (zeC) (2-1) 
d{z,T,k) 

where d{-,K) is the distance function from the closed subset K of the complex plane: 

diz, K) := inf Iz — t|. 

The function d{z,K) is non-negative, uniformly continuous on C so hk{z) is upper semi-continuous 

and /ifc(-z) = 00 on S^. 

Definition 2.1 A collection of background potentials 

VfciSfc ^ (-oo,cx)], A; = l,2,...,i? (2-2) 

is said to be admissible with respect to the (positive definite) interaction matrix A if the following 
conditions hold: 

[Al] the potentials Vk are lower semi-continuous on for all k, 

[A2] the sets {z e T,^ Vk{z) < 00} are of positive logarithmic capacity for all k, 

[A3] the functions 

^..(.,t):=Mf)±M)+,^,ln^ (2.3) 
are uniformly bounded from below, i.e. there exists an L G M such that 

Hjkiz,t)>L (2-4) 

on {{z, t) G Sj X Sfe : 27^ t} for all j, k = 1, . . . , R. Without loss of generality we can assume 
L = by adding a common constant to all the potentials so that 

Hjkiz,t)>0. (2-5) 

We will also assume (again, without loss of generality) that all the potentials are non- 
negative. 
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[A4] There exist constants < c < 1 and C such that (recall that a^k > 

H,k{z,t) > ^l^{V,{z) + Vk{t)) - ^. (2-6) 
The constant C can be chosen to be positive. 
[A5] The potentials are given such that the functions 

Qk{z):= (^yii^)+akMz)] = ^^Vk{z)+ ^kMz) (2-7) 

I: aki<0 ^ ^ «: a^iKO 

are bounded from below on (here < R — 1 is the number of negative a^i 's). 
Note that in the above sum k / / because of the assumption that a^f^ > 0. 

Definition 2.2 We define the weighted energy with interaction matrix A of a measure 
/2 = [/ii, . . . , /i/j] with Hj £ MiC^j) by 

R , R 



j,k k=l 

= Y jj Hjkiz,t)dfijiz)dfikit), (2-J 

j,k 

where Aii{K) stands for the set of all Borel probability measures supported on some set K C C 



Remark 2.1 The assumption [A|3] is a sufficient requirement to ensure that the definition of the 
functional ^(•) is well-posed and it is a rather mild assumption on the growth of the potentials 
near the overlap regions and infinity. Indeed ( with L = 0) 

^A,v(P) = Y.jj Hjk{z,t)disj{z)df,kit) > 0. (2-9) 

j,k 



Note also that if a conductor T,j is unbounded the condition (2-6) implies that 

^Vj{z) > a,j In |z - tol - l^i(to) - ^ (2-10) 

and hence Vj grows at least like a logarithm. In l2Uj the usual requirement is the stronger one that 

Vj{z)/ ln\z\ oo as z ^ oo. 

Remark 2.2 [A[4] is a stronger requirement which will be used for proving tightness (and therefore 
relative compactness) of a certain subfamily of measures over which ^(•) is guaranteed to attain 
its minimum value. 



Remark 2.3 [A|5] is yet stronger and assumes that all potentials have a suitable logarithmic growth 
near the common boundaries with those condensers carrying an opposite charge. This condition 
could be relaxed in some settings. 
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3 Existence and uniqueness of the equilibrium measure 



In this section we prove the existence and uniqueness of the equilibrium measure for the vector 
potential problem described above. Before stating our main theorem, we recall that a family of 
measures JT on a metric space X is called tight if for all e > there exists a compact set K <Z X 
such that ^{X \ K) < e for all measures /i G JT. The following theorem is a standard result in 
probability theory: 

Theorem 3.1 (Prokhorov's theorem [21J ) Let {X,d) be a separable metric space and M.i{X) 
the set of all Borel probability measures on X. 

• If a subset T C M.i{X) is a tight family of measures, then T is relatively compact in M.\{X) 
in the topology of weak convergence. 

• Conversely, if there exists an equivalent complete metric do on X then every relatively compact 
subset T of M-x^X) is also a tight family. 

We will use the following little lemma: 

Lemma 3.1 Let F:X^ [0, cxd] be a non-negative lower semi-continuous function on the locally 
compact metric space X satisfying 

lim F{x) = cx), (3-1) 

i.e. for all H > there exists a compact set K C X such that F{x) > H for all x ^ X \ K. Then 
for all H > iniF the family 

Th ■■= & Mi{X) : J Fdfi<H^ (3-2) 
is a non-empty tight subset of Aii{X). 

Proof. F attains its minimum at some point xq £ X since F is lower semi-continuous and 
liuix^oo F{x) = oo and therefore the Dirac measure dxg belongs to J^h- To prove the tightness of 
J^H, let e > be given. Since F goes to infinity "at the boundary" of X there exists a compact set 
K C X such that F{x) > ^ for all a; G X \ K. If G we have 

f,iX\K)= [ di,<^[ Fdi,<^ [ Fdfi<^H='-<e. (3-3) 
Jx\K Jx\K ^-n Jx ^ 

Q.E.D. 
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Define 

Ut{z) := an I In dfii{t), (3-4) 

k=l '' ' 

which is the logarithmic potential (external terms and self-potential together) experienced by the A;th 
charge component in the presence of fl only. 

Theorem 3.2 (see [20], Thm. VIM. 1.4) With the admissibility assumptions [Al] - [A5] above the 
following statements hold: 

1. The extremal value 

V^,^:=mil^y{fl) (3-5) 
of the functional Ij^^yi') is finite and there exists a unique (vector) measure /2* such that 

2. The components of jl* have finite logarithmic energy and compact support. Moreover, the Vj 's 
and the logarithmic potentials are bounded on the support of fi^ for all k = 1, ... , R. 

3. For j = 1, . . . , R the effective potential 

ip,iz):=uf{z) + Vjiz) (3-6) 

is bounded from below by a constant Fj ('Robin's constant^ on Sj, with the equality holding 
a.e. on the support of jjLj. 

Remark 3.1 The content of Thm. |3.2| is probably neither completely new nor very surprising and 
the proof is a rather straightforward generalization: the main improvement over the most common 
literature is the fact that we allow the condensers to overlap even if the corresponding term in the 
interaction matrix is negative. The assumption on the potentials that they provide a screening effect 
so that the equilibrium measures will not have support on the overlap region. The theorem will be 
instrumental in the proof of Thm. 5.1 which is the main result of the paper. 

Proof of Theorem 3.2[ First of all, we have to prove that 

by showing that there exists a vector measure with finite weighted energy. To this end, let fJ be the 
i?-tuple of measures whose fcth component rjk is the equilibrium measure of the standard weighted 
energy problem (in the sense of [20]) with potential Vk{z)/akk on the conductor for all k. (The 
potential Vk{z)/akk is admissible in the standard sense on since 

1 R 1 C 

Vk{z) - In \z\ > - In |z - to\ Vk{to) - In \z\ ^ oo (3-8) 

ttkk c ttkk cakk^ 
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as |2;| ^ oo for z G if is unbounded.) We know that % is supported on a compact set of the 
form 

^GS, :^<k4 (3-9) 



for some K^. G M. These sets are mutually disjoint by the growth condition (2-7) imposed on the 
potentials. The sum of the "diagonal" terms and the potential terms in the energy functional are 
finite for fj since this is just a linear combination of the individual weighted energies of the equilibrium 
measures rjk- The "off-diagonal" terms with positive interaction coefficient aui are bounded from 
above because the supports of r/^ and rji are separated by a positive distance; the terms with 
negative interaction coefficent are also bounded from above since rj^ and rn are compactly supported. 
Therefore 



Integrating the inequalities (2-6) it follows that 



R R 



^Ay(f^) = J2 [[ Hjk{z,t)disj{z)diskit) >{l-c)J2 [ Vk{z)disk{z) - a (3-11) 
j,fc=i -^-^ k=i 

We then study the minimization problem over the following set of probability measures: 

Y.J Vk{z)dfik{z) < (f3^C^Ay + '^ + ^)| ^-^1(^1) X ... x7Wi(Sfi) . (3-12) 

The extremal measure(s) are all contained in T since for a vector measure A ^ ^ we have 

I^yW > (1 - c) ^ / Vk{z)dXk{z) -C> V^^^ + 1. (3-13) 
k=i 

The function Yl,k^k{z) is non-negative, lower semi-continuous and goes to infinity as \z\ — > 00, and 
moreover 

^ (V^^ + C + 1)>0, (3-14) 



(1-c) 



hence, by Lemma 3.1 all projections of T to the individual factors is a non-empty tight family 



of measures. Using Prokhorov's Theorem 3.1 we know that there exists a measure /T* minimizing 
Ij^y{-) such that Jl^k=l^^* ^ ^- '^^^ existence of the (vector) equilibrium measure is therefore 
established. 

Note that now statement (2) follows immediately: indeed from the condition [S] that Hj i^ > 
(and also Vj > 0) it follows that 

\^^\dfi*,{z)dfit{t) + I 



^AV = «ii // lnr:-^df,*,{z)dfi\{t) + - I Vi{z)di^{{z) 
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> an jj In j^^dfi*{z)dfi\{t). (3-15) 

Thus the logarithmic energy of fi^ is bounded above by V^^/au- Repeating the argument for all 
/i*'s we have that all the logarithmic energies of the fi*'s are bounded above. 

On the other hand, these log-energies are also bounded below using (2-6) with j = k: 

II In j^^dfi*{z)df,*{t) > -| / V,{z)dfi*{z) - ^ (3-16) 

(boundedness from below follows since f Vj{z)d^j{z) is bounded above and appears with a negative 
coefficient in the formula). 

Now, using the fact that the quantities Hjk{z,t) are nonnegative due to (2-5) and condition 



( |3-12p it follows that 

cpjiz) = Vj{z) + ajk I In T^Mt) (3-17) 

is finite wherever Vj{z) is. Using condition [A[5] it also follows that it is lower semicontinuous, 
bounded from below on and hence admissible in the usual sense of minimizations of single 
measures [2D]. We also claim that ipj grows to infinity near all the contacts between and any 
for which ajk < 0. Suppose zq £ n (with ajk < 0); then on a compact neighborhood K of zq 
we have 

^j{z) > Vj{z) + ajkhkiz) + Mk (3-18) 



for some finite constant Mk (which -of course- depends on K). From ^ then 

Vjiz) + ajkhkiz) + Mk> ^^^^^A^) + (3-19) 

k^j 

where sj < R \s the number of negative ajk {j 7^ k). Since Vj{z) tends to infinity at the contact 
points (from the same condition [A|5]) then so must be for (pj. 
Note also that 

j 

and hence (as in [20]) each single nj^i, is the minimizer of the single variational problem on Ej under 
the effective potential ipj. From the standard results it follows that the support of jj.* is contained in 
the set where ipj is bounded, which, due to our assumptions, are all compact and at finite nonzero 
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distance from the common overlaps. This proves that the components of jl* are actually compactly 
supported. 

Uniqueness as well as the remaining properties are established essentially in the same way as in 
[20J, Thm. 1 Chap. VIII using the positive definiteness of the interaction matrix^, which guarantees 
the convexity of the functional. 

Q.E.D. 



4 The special case 



We now specialize the above setting to the following collection of R conductors: 
and interaction matrix 




-52^3 



.- (-1)^-1 [0,oo) (j = l,2,...,i?), 
2g? -qiq2 



(4-1) 



A :-- 



-qm "^^ql -q2q3 



















(4-2) 



Under the assumptions on the growth of the potentials Vj{x) near the only common boundary 
point X = 0, Thm. |3.2| guarantees the existence of a unique vector minimizer. 



We now investigate the regularity properties under the rather comfortable assumption that the 
potentials Vj are real analytic on \ {0} for all j; this is in addition to the host of assumptions 



specified in Def. 2.1 



In order to simplify slightly some algebraic manipulations to come we re-define the problem by 
rescaling the component of the vector of probability measures fij i— > qjfij so that now the interaction 
matrix becomes the simpler 



(4-3) 





2 


-1 


. 


.. 




-1 


2 


-1 . 


.. 


A := 





-1 


2 . 


. . 










. 


. . 2 



The electrostatic energy can be rewritten as 

^A,v(P) = 2 J]] // In— — (i/ij(x)(i^j(y) - ^ 



F - y\ 



Ft « 

+2^ / Vj{x)dfxj{x). 



In, — ^ — -dfj.j{x)dfij+i{y) (4-4) 
\x-y\ 



(4-5) 
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As explained in the previous section, the above minimization problem has the interesting property 
that the same equilibrium measure is achieved by minimizing only one component of it in the mean 
field of the neighbors and, moreover, the supports of the minimizers satisfy 

supp(pj) n supp(pj+i) = . (4-6) 



Corollary 4.1 Let fl be the vector equilibrium measure for the above problem. For any 1 < k < R 
we have that 

/ ln^-d/i(z)dMt)+2 /" Vk{z)dfi{z) (4-7) 
is minimized by the same fik> where the effective potentials Vk are 

Vi{z) := \vi{z)-\l \ny^d^i2{t) (4-8) 

■■= \Vk{z)-\f \ny^diJik+i{t)-\ f \n-^diJik_i{t) (4-9) 

Vr{z) := \vR{z)-\j \ny^d^iR.^{t). (4-10) 

Note that the effective potential differs from the original potential by harmonic potentials because 
the supports of iXk±i are disjoint from the support of //fc. 
We recall the following theorem: 

Theorem 4.1 (Thm. 1.34 in [16| ) If the external potential belongs to the class , k > 3 then 
the equilibrium measure is absolutely continuous and its density is Holder continuous of order ^. 



Combining Cor. |4.1| with Thm. |4.1| we have that the solution of our equilibrium problem consists 
of equilibrium measures which are absolutely continuous with respect to the Lebesgue measure with 
densities pj at least Holder-^ continuous as long as the external potentials are at least C^. Moreover 
the supports of these equilibrium measures have a finite positive distance from the origin. 

Our next goal is to prove that the supports of the pj's consist of a finite union of disjoint compact 
intervals. For that we need a pseudo-algebraic curve given in the next section. 

5 Spectral curve 

Since the equilibrium measures have a smooth density we can now proceed with some manipulations 
using the variational equations. 

For the remainder of the paper we will make the following additional assumption (besides those 



n Def. 2.1) on the nature of the potentials Vj: 
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Assumption: the derivative of tlie potential V- is the restriction to S° := (— 1)-^~^(0, cxd) of 
a real analytic function defined in a neighborhood of the real axis possessing at most isolated polar 
singularities on R \ S j . 

For a function / analytic on C \ T, where F is an oriented smooth curve, we denote 

S{f){x):=U{x)+Ux) , A(/)(x) :=/+(x)-^(x) , xGT. (5-1) 

where the subscripts denote the boundary values. We remind the reader that under our assumptions, 
the equilibrium measures satisfy eq. (4-6). 

Definition 5.1 For the solution p of the variational problem, we define the resolvents as the 



expressions 



Wj{z) :-- 



Pj{x)dx 

Z-X 



z £ C \ supp(pj 



(5-2) 



The variational equations imply the following identities for j = 1, . . . , R: 

S{Wj){x) = V;{x) + Wj+i + Wj.i 
A{Wj){x) = —2i-Kpj{x), X £ supp(pj) 



(5-3) 



where we have convened that Wq = TV/j+i = 0. Note that, under our assumptions for the growth 
of the potentials Vj at the contact points between conductors (in this case the origin), the support 
of pj is disjoint from the supports of pj±i and hence the resolvents on the rhs of the above equation 
are continuous on supp(pj). 

The following manipulations are purely algebraic: we first introduce the new vector of functions 



Yi 



Yr 



-1 



1 









' Wi ' 








+ 




j 


< 











(5-4) 



(-1)« 

Trivial linear algebra implies then the following relations for the newly defined functions Yj: 



S{Y) 
S{Y2) 
S{Ys) 

S{Yn-i) 
S{Yr) 



-Y2 

-Yi - Ys 

-Y2 - n 



A(n) 

A(l2) 
A(l3) 



2i7rpi 
-2i-Kp2 

2iTTp3 



-Yr_2-Yr A(Yn-i) 
-Yn-i A{Yr) 



(-l)«2i7rp/j_i 
{-l)R+^2i7rpn 



on supp(/9i) 
on supp(/92) 
on supp(/33) 

on supp(/9K-i) 
on supp(pi?). 



(5-5) 



The above relation should be understood at all points that do not coincide with some of the isolated 
singularities of some potential Vj (points of which type there are only finitely many within any 
compact set). 
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Define then the functions 



Zo:=yi, Zi:=-Yi-Y2, Z2 := Fs + ils , • • • , ^R-i = , Zr := {-1)''Yr. 

(5-6) 

Then 

Proposition 5.1 All symmetric polynomials of {Zj}Q<j<R are real analytic in the common domain 
of analyticity of the potentials, namely they have no discontinuities on the supports of the measures 

Pi- 

Proof. A direct algebraic computation using the boundary values of the {Yj} functions gives the 
following boundary values of the functions Zj-. 

2Zo± = -Y2 ± 2m pi (5-7) 

= i -y2^2i-Kpi = 2Zo=p on supp(pi) 

\ -Yi +Y^± 2iTTp2 on supp(p2) " 

^ f -Yi+Y3=f 2i7rp2 = 2Zi^ on supp(p2) ,g_gs 
\ Y2-Y4± 2i7rp3 on supp(p3) 



(5-10) 



-1)^ 1(-Yr_2 + Yr) =f 2iTTpR_i = 2Z(^_2)^ on supp{pR^i] 



-\{-lr-'YR_,~±2^.pR " onsupp(p^) ('"''^ 

2Zr^ ={-lf-^YR^i^2iTrpR = 2Z^R_i)^ on supp(pij) (5-12) 

Consider a symmetric polynomial Pk ■= 2^ {Zq^ -|- . . . -|- Zr^^ and its boundary values on, say, 
supp(/9i); we see above that Zq^ = Z\^ and hence Z^ -|- Z^ has no jump there. The support of 
P2 has no intersection with Si and supp(pi) (see (4-6)) due to our assumptions, and hence Z2 may 
have a jump on supp(/3i) only if the support of pz has some intersection with it. In that case anyway 
Z2^ = Z'i^ and hence also Z2 -|- Z^ has no jump on supp(p3) n supp(/9i). 

In general on supp(p,fc) nsupp(pi) we have Z^^ = and so the same argument apply. In short 
one can thus check that all the jumps that may a priori occur in fact cancel out in a similar way. 

Repeating the argument for all the other supp(pj) instead of supp(pi) proves that the expression 
has no jump on any of the supports, and since a priori it can have jumps only there, then it has no 
jumps at all. Invoking Morera's theorem, we see that the symmetric polynomials of the Z^s can be 
extended analytically across the supports of the pj's. 

Finally, the statement that the symmetric polynomials are real analytic follows from the following 
reasoning: the Zj's are linear expressions in the WjS and the potentials. In particular they are 
analytic off the real axis (where all the WjS are) and in the common domain of analyticity of the 
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potentials. The same then applies to the symmetric polynomials in the Zj's. Finally, on an open 
interval in M, as long as it is outside of all the supports of the vector measure, the Zj are all real 
analytic functions since Wj's are. This concludes the proof. 

Q.E.D. 

A consequence of this proposition is that 
Theorem 5.1 The functions are solution of a pseudo-algebraic equation of the form 

z^+^ + C2{x)z^-^ + ... + CR+lix) = (5-13) 

where Cj{x) := (—I)-' J^e^ e- ' ' ' '^ij ^''^ (real) analytic functions on the common domain of 
analyticity of the potentials. 

Remark 5.1 This result is the direct analogue of the results about the existence of the spectral 
curve for the one-matrix model 117] which was established on a rigorous ground in fT^ . In a different 



context of matrix models with external source Thm. 5.1 is conceptually similar to the result in /7^ . 
Proof of Thm. Hi] We set 

R 

E{z,x):=ll{z-Zj{x)) , (5-14) 

j=0 

and expand the polynomial in z. Clearly we have Zq + Zi + . . . + Zr = and hence the coefficient 
Ci vanishes identically. The other coefficients are polynomials in the elementary symmetric functions 
already shown to be real analytic and hence sharing the same property. 

Q.E.D. 



Corollary 5.1 The densities pj are supported on a finite union of compact intervals. Moreover 

the supports of pj and pj±i are disjoint. 

Proof. 

The supports of the measures are in corre- 
spondence with the jumps of the algebraic so- 
lutions of E{z,x) = 0; in particular the set of 
endpoints of the supports must be a subset of 
the zeroes or poles of the discriminant that be- 
long to M. Since the only singularities that these 
may have come from those of the derivatives of 
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Figure 1: The Hurwitz diagram of the spectral curve. 



the potentials Vj{x) on the real axis, and these 
have been assumed to be meromorphic on M and 
be otherwise real analytic, then the discriminant 
of the pseudo-algebraic equation cannot have 
infinitely many zeroes on a compact set. We 
also know that the measures pj are compactly 
supported a priori and hence there can be only 
finitely many intervals of support. 



Q.E.D. 



Putting together Prop. 5.1 and Thm. 5.1 we can rephrase the properties of the functions Zj{x) 



by saying that they are the R+1 branches of the polynomial equation (5-13), thus defining an 
{R + l)-fold covering of (a neighborhood of) the real axis. The neighborhood is the maximal 
common neighborhood of joint analyticity of the potentials Vj{x). The various s/7eets defined by the 
functions Zj{z) are glued toghether along the supports of the equilibrium measures pj in a "chain" 
of sheets as the Hurwitz diagram in Fig. [T] shows. 

Remark 5.2 In Jdj a similar problem was considered in the context of multiple orthogonality for 
Nikishin systems on conductors without intersection and with fixed weights: this corresponds to 
the case of a minimization problem without external fields. It was shown that an algebraic curve 
similarly arises; in the formulation of [6] the algebraic curve involves, rather than the resolvents, their 
exponentiated antiderivatives ^'j 's, namely 

Wj = ^In^ J x) (5-15) 
dx 

and a mixture of algebraic geometry and geometric function theory was used to investigate their 
properties. In particular the functions figured in an algebraic equation (see eq. (2.1) in f^) as 
the various determinations of a polynomial relation 

^^^^ +r^{x)^^ + . . . + r^{x)^ + r^^^{x) =0 , rj£C[x] (5-16) 

with the discriminant (w.r.t. "i>) vanishing at the endpoints of the supports for the measures of the 
corresponding Nikishin problem. Along similar lines, examples of curves of algebraic type for Nikishin 
systems with special choices of external fields were recently obtained in 
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6 An explicit example 



We consider the case with R = 2 and the two potentials are the same Vi{x) = V2{—x) and are of 
the simplest possible form that satisfies our requirements 

Vi{x) = bx — alnx , x>0; V2{x) = —bx — aln{—x) , x < (6-1) 

where both a,b > 0. 

Quite clearly we can rescale the axis and set 6 = 1 without loss of generality. 



Using the expressions for the coefficients of the spectral curve (Thm. 5-13) in terms of the 
potentials V\ = V and V2 = V* = V{—x) we have 

E{z,x) = - R{x)z - D{x) = (6-2) 

where, on account of the fact that the derivative of the potentials have a simple pole at x = 0, the 
coefficients R{x),D{x) have at most a double pole there. From the relationship between the three 



branches of Z and the resolvents Wi, W2 (eq. 5-4) we have 



Z^^){x) = -Wi-- + \ (6-3) 
X 6 

Z^^\x) = W2 + - + \ (6-4) 
X 6 

Z(i)(x) = -ZW(x) - Z(2)(x) = Wi{x) - W2{x) + — (6-5) 

X 

and hence the general forms that we can expect for the coefficients of the algebraic curve are 

R{x) = + + - 
3 X 

where the constants A,B,C have yet to be determined. 

The spectral curve — Rz — D = has in general 5 finite branchpoints (which is incompatible 
with the requirements of compactness of the support of the measures) and requiring that there are 
< 4 branchpoints and symmetrically placed around the origin(by looking at the discriminant of the 
equation) imposes that B = C = 0. 

The ensuing spectral curve is 

, /I a^x f2a'^ + 3A 2\ , 
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and a suitable rational uniformization of this curve is 



3a2 (A2-(l + ^/a2))a4 

The three points above a; = oo are A = ±1, and Z is regular there. 

We see that the condition that the measures pi,p2 have unit mass requires that 

res Z^^^dx = 1 + a , res Z^^^dx = -1 - a . (6-10) 

x=co x=co 

We need only to decide which point A = ±1,0 correspond to the three points over infinity. But this 
is achieved by inspection of the behavior of Y{\) and -'^(A) near the three points A = 0, 1, —1. 0. 
By this inspection we have 

A = 1 ^ ooi (6-11) 
A = -1 ^ 002 (6-12) 
A = ^ ooo . (6-13) 

Computing the residues of Zdx = ZX'dX at these points we have 

res Z(°)da; = Va'^ + A = 1 + a (6-14) 

x=co 

res Z^'^^dx = -Va^ + A = -1 - a (6-15) 

x=oo 

which imply that ^ = 2a -|- 1. 

Collecting the above, we have found that 

a + 1 2a + l ( 1 1 \ 

^ = — -2^1aTT+a^J ^'-''^ 

2a2 + 6a + 3 (2a + l)(a+l) 

3a2 (A2-((a+l)Va2))a4 ^ > 

and the algebraic equation for z = Z{X) in terms of x = X{\) becomes 



3 A , a2\ /2a2 + 6a + 6 2 \ ^ ,^ 

It is possible to write explicitly the expressions of the branchpoints in terms of a but it is not very 
interesting per se, except to discuss their behaviors in different regimes of a; we find that for a > 
there are four symmetric branchpoints on the real axis and the inmost ones tend to zero as a — > 0, 
whereas they all tend to infinity as a — >^ oo according to ib(a it 2-y/a) -|- 0{1). 
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Figure 2: Some examples for the equilibrium measure for the example worked out in the text, and 
a = 0, 1, 2, 3 respectively from left to right. In red is the graph of the potential Vi. The symmetry implies 
that the other equilibrium measure is simply the reflection of this around the ordinate axis. The units 
for the axes are the same in all cases. The growth of the density at x = for a = is 0(a;~^/'^). Near 
the other edges the vanishing is of the form 0{{x — a)^). 



It is interesting to note that for a = our 
general theorem does not apply: the potentials 
are finite on the common boundary of the con- 
densers and hence cannot prevent accumulation 
of charge there. However the algebraic solution 
we have obtained is perfectly well-defined for 
a = giving the algebraic relation 

o z 2 2 

z-^ 

3 
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,^2 



A short exercise using Cardano's formulae shows that the origin is a branchpoint of order 3 and thus 
corresponding to the Hurwitz diagram on the side. 

The behavior of the equilibrium densities pj near the origin is (expectedly) x~3. 



7 Concluding remarks 

We point out a few shortcomings and interesting open questions about the above problem. 

The first problem would be to relax the growth condition of the potentials near common points 
of boundaries, if not in the general case at least in the specific example given in the second half of 
the paper, where we consider conductors being subsets of the real axis. 

The importance of this setup is in relation to the asymptotic analysis of certain biorthogonal 
polynomials studied elsewhere [8J and their relationship with a random multi-matrix model [lOj. 

In that setting, having bounded potentials near the origin G M would allow the occurrence of 
new universality classes where new parametrices for the corresponding 3 x 3 (in the simplest case) 
Riemann-Hilbert problem would have to be constructed. 
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Based on heuristic considerations involving the analysis of the spectral curve of said RH problems, 
the density of eigenvalues should have a behavior of type x~3 near the origin (to be compared with 
,T~2 for the usual hard-edge in the Hermitian matrix model). Generalization involving chain matrix 
model would allow arbitrary —2 behavior, p < q. However, for all these analyses to take place the 
corresponding equilibrium problem should be analyzed from the point of view of potential theory, 
allowing bounded potentials near the point of contact. 
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